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Abstract 

We study static cylindrically symmetric vacuum solutions in Weyl coordinates in the context of 



O 
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(^T) ■ the metric f(R) theories of gravity. The set of the modified Einstein equations is reduced to a single 



equation and it is shown how one can construct exact solutions corresponding to different f(R) 
models. In particular the family of solutions with constant Ricci scalar (R = Rg) is found explicitly 
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iJ^j" which, as a special case (R = 0), includes the exterior spacetime of a cosmic string. Another new 

■ solution with constant, non-zero Ricci scalar is obtained and its possible relation to the Linet-Tian 

> . 

C*~) , solution in general relativity is discussed. 
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I. INTRODUCTION 

As pointed out nicely by Weinberg in his seminal paper [l| on the cosmological constant 
problem, " Physics thrives on crisis" . Perhaps the biggest crisis of the 20-th century physics 
which was carried over to the 21-st century is the so called cosmological constant or dark 
energy problem: the 120 orders of magnitude difference between the observational and 
theoretical values of the vacuum energy density. To overcome this crisis, different proposals 
have been put forward with the hope of obtaining a consistent theoretical background to 
the recent observation of an expanding universe which is seemingly not accessible through 
the standard model of cosmology. Obviously one might trace back this lack of a proper 
theoretical explanation to the basic ingredients of the standard model of cosmology, one of 
which being the Einstein field equations derived from the Einstein-Hilbert action. Modified 
or alternative theories of gravity is the paradigm under which all those theories which differ 
from the Einsteinian gravity are studied. One of these non-Einsteinian theories, that arose 
a lot of enthusiasm recently, is the so called f(R) gravity in which a function, f(R), replaces 
the Einstein-Hilbert (gravitational) Lagrangian R. It seems that f(R) actions were first 
contemplated by Eddington jsj and later on rigorously studied by Buchdahl in the context 
of nonsingular oscillating cosmologies. These theories could be thought of as a special kind 
of higher derivative gravitational theories. Having shown that these models are equivalent 
to scalar-tensor models of gravity it is obvious that one should in the first place check their 
consistency with the solar system tests of Einstein gravity 4j. In most of the models it is not 
possible to satisfy these tests and at the same time to account for the accelerated expansion 
of the Universe without bringing in new degrees of freedom. Recently it was shown that there 
are models of f(R) gravity in which one could account for both the accelerated expansion of 
the universe and the solar system tests without introducing the cosmological constant 
This is why spherically symmetric solutions are the most widely studied exact solutions in 
the context of f(R) gravity. Apart from this it is interesting, at least from theoretical point of 
view, to consider other exact solutions of the modified Einstein equations of f(R) theory. As 
in the case of the Einstein-Hilbert action one could derive field equations in f(R) gravity in 
two different approaches, the so called metric and Palatini approaches. But in f(R) actions, 
unlike the Einstein-Hilbert action or its modified version (one with a cosmological constant 
term), the field equations obtained by the two approaches are not the same in general. In 
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what follows we will be interested only in metric f(R) theories of gravity in which connection 
is dependent on metric g^ v with respect to which, the action is varied. In ordinary GR there 
are not that many exact solutions of the field equations for a given symmetry. Being higher 
derivative theory it is not unexpected to find more exact solutions in f(R) gravity and this 
turns out to be the case for spherically symmetric solutions 7|. Since cylindrical symmetry 
is the next symmetry considered normally in the study of exact solutions in GR (not just for 
theoretical reasons but also because they might have physical realization in objects such as 
cosmic strings) it seems natural to extend the studies of exact solutions in f(R) theories in the 
same direction. Looking for solutions with a different symmetry, as a first step, we consider 
static cylindrically symmetric vacuum solutions of the f(R) modified Einstein equations in 
this letter. It is shown how one can reduce the set of equations into a single equation which 
could then be utilized to construct explicit solutions. For constant curvature solutions, using 
the general form of a cylindrically symmetric solution in Weyl coordinates we find, among 
possible solutions, a generalized form of a conical (zero curvature) spacetime as well as two 
new (non-zero curvature) solutions with one of their parameters chosen so that it is related 
naively to the cosmological constant (section IV). We discuss possible relation to the so 
called Linet-Tian (LT) solution of the modified Einstein field equations which is the 

cylindrical analogue of the Schwarzschild-de Sitter spacetime. 

II. FIELD EQUATIONS IN /(R) GRAVITY 

In this section we give a brief review of the field equations in f(R) gravity. The action 
for f(R) gravity is given by 

S = J{f{R) + Z m )^-gd A x. (1) 

The field equation resulting from this action in the metric approach, i.e. assuming that the 
connection is the Levi-Civita connection and the variation is done with respect to the metric 
SV, is given by 

1 _ g T l 

where the gravitational stress-energy tensor is 



— R^iv n R9fj,u — T^j, + T _ I/ D ^ , (2) 



with F(R) = df(R)/dR and T™, the standard matter stress-energy tensor derived from 
the matter Lagrangian £ m in the action (1). The vacuum equations of motion, i.e. in the 
absence of matter, are given by, 

F(R)R flu - l -f{R)g^ - V,V U F(R) + g^UF{R) = (4) 

Contraction of the above field equations gives the following relation between f(R) and its 
derivative F(R) 

F(R)R — 2f(R) + 3nF(R) = 0, (5) 

which will be employed later both to simplify the field equations and to find the general 
form of the f{R) function. 



III. CYLINDRICALLY SYMMETRIC VACUUM SOLUTIONS 

Interested in the static cylindrically symmetric solutions of the vacuum field equations 
(4), we start with the general form of such a metric in the cylindrical Weyl coordinates 
(t,r,<p,z) given by [lo|; 

9lw = diag{-e 2k ^~ 2 < r \ e 2k ^~ 2 < r \ W ( r f e -Mr)^ e 2 «M). (6) 

The corresponding scalar curvature is 

„ —wu" + wk" — w'v! + w" + wu' 2 , , 

R= > 7 

we 

in which ' = j-- Using equation (5), the modified Einstein equations become 

FR„ V - V,V U F = l -g, v {FR - UF{R)). (8) 

As in the spherical case, since the metric only depends on the cylindrical radial coordinate r, 
one can view equation (8) as a set of differential equations for functions F(r), u(r), k(r) and 
w(r). In this case both sides are diagonal and hence we have four equations. Differentiating 
equation (5) with respect to r we have the extra, consistency relation for F(r), 

RF' - R'F + 3(DF)' = 0. (9) 
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Any solution of equation (8) must satisfy this relation in order to be also a solution of the 
original modified Einstein's equations. From equation (8) it is found that 

FR^ - V M V M F = h FR _ UF{R)). (10) 
9fj.fi ^ 

In other words the combination = F ' Rw '~ Va ' Va " f (with fix indices) is independent of the 
index fi and therefore = A u for all fi, v. This allows us to write the following independent 
field equations; 

_ p" + 2F'(k' -u')+F(-— + ^ + 2u 2 ) =0 (11) 

w w 

J I I II 

Fw 2 (-k + — +F fcw 2 -w =0 12 

w w 

„ , ,n ii kw 2w u' . i, i i . . 

Fw(-k +2u + )+F (wk - 2uw) = 0, 13 

w w 

corresponding to A t = A r , A t = A^ and A t = A z respectively. Therefore, any set of functions 
u(r), k(r) and and w(r) satisfying the above equations would be a solution of the modified 
Einstein field equations (8) for a given F(r) satisfying equation (5). Obviously it is not an 
easy task to find a general solution to the above equations, so in the following section we 
discuss the simple but important case of solutions with constant curvature. 

IV. CONSTANT CURVATURE SOLUTIONS 

It is known that some of the vacuum constant curvature solutions in f(R) gravity are 
equivalent to vacuum solutions in Einstein theory with the same symmetry. For example 
it is shown in [?| that in the spherically symmetric case the corresponding f(R) solutions 
include the Schwarzschild-de-Sitter space for a specific choice of one of the constants of 
integration. For cylindrical symmetry, in Einstein gravity, static vacuum solutions were 



found almost immediately after their spherical counter parts by Levi-Civita ll| but those 
with a cosmological constant have to wait another 60 years to be found independently by 
Linet 8j and Tian 9|. Their solution reduces to that of a cosmic string in the limit r — > 0. 
Looking for cylindrically symmetric solutions in f(R) gravity, here we consider the simple 
but physically important case of static constant curvature spacetimes. To do so, taking 
R = constant in the field equations (11), (12) and (13) , we arrive at the following set of 
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equations: 



2u" + ^- 


n k'w' 


= 


w 


w 




k" + kv 


w" 


= 


w 


w 




^w'u' k'w' 


- k" 


= 0. 




w w 







(14) 
(15) 
(16) 



From equations (15) and (16) one could obtain the following two equations: 



u'= l -^l (17) 

k> = ^i, (18) 
w 

in which ci and c 2 are integration constants. By substituting equations (17) and (18) in (14) 
we obtain the following differential equation for w(r): 

1 w' + c 2 2 w" _ 2 w'(d + w') 

2 w w " w 2 

One could solve the above equation by inspection and the solutions discussed below corre- 
spond to solutions with zero and non-zero values of the Ricci scalar R. 

Case (1) : solution with R = 

Looking at equation (19) one could obviously arrange for a solution of w(r) linear in r 
(i.e w" = 0) in which case the metric functions are given as follows; 

m = c 3 ± a /— \n(w) (20) 

V c 6 

fc = C4 + -ln(-) (21) 
w = c 6 r + c 7 . (22) 

It could be seen that this is a Ricci flat solution (i.e R — in (7)) in which we should identify 
which one of the constants C3 to C7 correspond to physical parameters of the spacetime and 



161 ] . Substituting the above 



which ones could be absorbed into the coordinate redefinitions 
functions back into the metric form ([6]) we obtain 

-2(c 3 ± /|pnp) ( 2{c i+ ^lnp),dp 2 2 n 2 \ 2(cj±jf fcv) , 2 



ds 2 = e~^V ( e ^~^(p_ _ dt ^ + p* d(f ? ) + ^'dz 2 , (23) 
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in which C4 = C4 — ^Inc^ and p = w = c$r is the new radial coordinate by setting C7 = 
without any loss of generality. One can also show that through the following redefinitions 
of the constants and the coordinates 



™=y| (24, 

t = e £4 - C3 A-(4i)+it (25) 

z = e c 3 Am(m^i)+i z (26) 

(j) = e~ C3 A rnimkn+l (j) (27) 

P = ^m(m T l) + lp (28) 
e 4~ c 3 

A = e c e , (29) 



the above metric reduces to 



and on applying the following complex transformation [h]] 

t — > iz z — > it, (31) 
it transforms into the following well known metric 

ds 2 = p^TD ( df) 2 +d f) + p^2m d ^ _ p±2m d f j (33) 

which is formally similar to the Levi-Civita's static cylindrically symmetric solution in GR 
normally written without ± sign but with the constant m taking both positive and negative 
values. It should also be noted that the range of the variable <ft is not in general (0, 2tt], not 
even for the flat cases of m = 0, 1. 

In the case of m = the spacetime is conical with a deficit angle corresponding to the 
exterior metric of a cosmic string with the following line element [12] 

ds 2 = (dp 2 + dz 2 ) + a 2 p 2 d(j) 2 - dt \ (33) 

in which oq = CQe~ Ci is the conical parameter related to the gravitational mass per unit 
length of the spacetime, r\, as [3] 

a = c 6 e- C4 = 1-477, (34) 



such that < ao < 1 for < C4 < oo (taking cq = 1). This metric, exposing the geometry 
around a straight cosmic string, is locally identical to that of flat spacetime however it 
is not globally Euclidian since the angle = ao0 varies in the range < < B where 
B = 2Tta = 2tt(1 — 477). 

From equation (5) it is seen that this metric is a solution for any form of f(R) for which 
f(R = 0) = 0, in other words, f(R) must be a linear superposition of R n with n > 1 of 
which the commonly considered model is f(R) = R + ^t- and obviously it is not a solution 
of the widely studied model for which f(R) = R — 

Case (2) : solutions with R = constant 7^ 

We discuss two solutions of this type here: 
A)-The first solution could be obtained by noting the simplifying fact that equation (16) is 
satisfied for k = 2u so that equations (14) and (15) reduce to, 

2u- + ^ - 4^ = (35) 
w w 

n „ 2u'w' w" n 

2u" + = 0, (36) 

w w 

leading to the following equation, 

/ + u '2_^ =0 , (37) 
w 

and consequently to the following relation between functions u and w, 

e u = D J wdp, (38) 

in which D is a constant to be determined later. On the other hand using (I3T|) . from equation 
(7) we have, 

w" 

R = R = —e~ 4u . (39) 
w 

Substituting from fl38|) we have; 

— = Ro(D [ wdp) 4 (40) 
w J 

. A solution to this integro-differential equation is given by; 

-"" 3/2 D ^tk )l " (41) 
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so that the metric takes the form 

ds 2 = AD 2 p- 1 (-dt 2 + dz 2 + dp 2 ) + f^d4> 2 , (42) 

or in a new coordinate system with p = ADp 1 / 2 

ds 2 = 6AD 4 p- 2 {-dt 2 + dz 2 ) + dp 2 + QAD 2 p- i d(f) 2 . (43) 

Now if one takes Rq = 4A, as in the case of the Einstein field equations in the presence of 
A, then the constant D introduced above is related to the cosmo logical constant through 

D 4 = ii A -i. (44) 
256 v ; 

Of course one should be careful with this kind of identification of the spacetime parameter 
as will be discussed later. 

B)-The other solution with constant non-zero Ricci scalar could be found by starting from 
the definition p := p(w) = ^ in terms of which the equation (19) can be integrated to 
obtain the following first integral: 

ln(w) — - \n(3p(w) 2 — 2p(w)c 2 — c 2 , + Ap(w)ci) 
-{c 2 \ arctan( 2 ^= 2 ° 2 + ° l ) + 

6 V~ C 2 + °2Cl ~ Cj 

c x \ axctan(j ^4^ ^ + C \ ))(J-c 2 + c 2Cl - c 2 )-' + C = 0. (45) 
6 V~ C 2 + c 2Ci - q v 

Also in terms of the same function the radial coordinate and the metric functions are given 
by; 

dw 



p(w) 



(46) 



u(w) = 1/2 f P -^l±^dw (47) 
J wp(w) 

k(w) = [ P -^±^dw. (48) 
J wp(w) 

Again it does not seem to be an easy task to find solutions for the complicated equation 
(H5|) . but in principle for each set of the values for the integration constants ci, c 2 and C we 
have a solution for p(w) and correspondingly solutions for metric functions k(r) and u(r). 
One such solution, looking at equation (JI5]) . could be obtained by choosing c 2 = 2c\. In this 
case, defining A = e c , we find the following solution for the metric functions 

1 f \/A i w' i + 4ci 2 \ 

k(w) = In (w) ■= arctanh I J (49) 

v3 V 2ci / 



w(w) = — In yw) y= arctanh I I , (50) 

in which w is a solution to the following equation 

3(^)2 _ e 3C w 3_ 4ci 2 = _ (51) 

dr 

1 /3 

Apart from the trivial solution u; = (— 4ci 2 ) e~ c (which is indeed equal to zero by equation 
(19) for a constant w and the fact that C2 = 1c\ ) one could show that there is a solution in 
terms of the Weierstrass p function [17| as follows; 

2 l/3 

iy = p^eierstrassp( — v / 3e c r + d, 0, -4ci 2 )e~ c 2 2/3 . (52) 

The constant curvature of the spacetimes given by the functions (j4"9l) . fl50l) and (1521) could be 
found through equation (7) in the coordinate system (t, w, <fr, z), for which we find Rq = A 3 . 
Now if again this is going to be compared with the solutions of the Einstein field equations 
in the presence of the cosmological constant for which Rq = 4A, then the correspondence 
Ro = e 3C = 4A will fix the value of the constant C as 

C = hn4A. (53) 

By looking at equation (5) it is clear that for constant Ricci scalar solutions (R = Rq), no 
matter what the symmetry, f(R) should satisfy the relation F(R )Rq = 2f(R ). Unlike 
the solution with R = 0, the above two solutions satisfy this relation for the commonly 
considered model of f(R) = R — ^ with R 2 = 3/i 4 if we have 



5 ,o 3 

1 



D 8 = (— ) 2 ^ (54) 
v 64 ; /i 4 



C = -/n(3/i 4 ) (55) 
6 



in cases A and B respectively 



V. DISCUSSION AND SUMMARY 



Studies on the exact f(R) gravity solutions are mostly restricted to the spherically sym- 
metric case mainly due to the solvability of the equations and also more importantly the fact 
that one could compare the results with the solar system observations/experiments based 
on the schwarzschild solution as the spacetime metric around a spherical mass such as Sun. 
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Here we have examined cylindrically symmetric solutions in metric f(R) gravity in a general 
form given by (6). Restricted by the complexity of the field equations we have only exam- 
ined solutions with constant Ricci scalar. In the case of Ricci flat solutions we have found 
a generalized form of a conical spacetime (with zero curvature) which, as a special case, 
includes the cosmic string spacetime. In the non-zero Ricci scalar case we have obtained 
two new solutions and in both of them a parameter is identified as the cosmological constant 
through the comparison of their Ricci scalar with that of the modified Einstein field equa- 
tions (in the presence of the cosmological constant). Obviously neither of these spacetimes 
are asymptotically flat nor they behave regularly in the limit p — > 0, and so one can not 
use their asymptotic behaviour to compare their parameters with those solutions known in 
GR having the same symmetry. Therefore one should note that either of the identifications 
(jHJ) and fl53l) , is a very naive one in the sense that we have not compared the corresponding 
spacetimes (nor we have studied their specifications) with one already known in the context 
of Einstein GR. Even if there are any solutions in GR comparable to these solutions (which 
we are not aware of), since the field equations in f(R) gravity are in general of higher order 
compared to their counterparts in GR, the correct identification should be made through 
matchings of different patches of the whole manifold. For example in the cylindrical case 
these might correspond to the exterior and interior solutions of a cylindrical shell if one of 
the parameters is going to be interpreted one way or another as the mass of the shell [?J ■ 

re metric fl43l) 



Another point need to be mentioned Incidentally it should be noted that t 
has the same general form as the LT solution introduced in Q (see also Q) [IS]. Unless 
one could find a solution to the equation (1401) which exactly corresponds to the LT solution 
in Einstein gravity, it seems a reasonable conjecture to say that the solution (|43j) is the 
f(R) analogue of LT solution in GR. Finally it should be noted that in the vacuum case, 
or R = constant the metric and the Palatini f(R) formulations are dynamically equivalent 



3- 



There 



same f(R) (l9| 



ore all the above solutions are also solutions in the Palatini gravity with the 
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